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The spherical-to-prolate shape transition in neutron-rich Cr isotopes from A'^ = 34 to 42 is stud- 
ied by solving the collective Schrodinger equation for the five-dimensional quadrupole collective 
Hamiltonian. The collective potential and inertial functions are microscopically derived with use 
of the constrained Hartree-Fock-Bogoliubov plus local quasiparticle random-phase approximation 
method. Nature of the quadrupole collectivity of low-lying states is discussed by evaluating exci- 
tation spectra and electric quadrupole moments and transition strengths. The result of calculation 
indicates that Cr isotopes around ^''Cr are prolately deformed but still possess transitional character; 
large-amplitude shape fluctuations dominate in their low-lying states. 
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I. INTRODUCTION 

Recent experiments on neutron-rich Cr isotopes show 
that quadrupole collectivity appreciably develops toward 
64Cr with iV = 40 [iHi. Going from ^^Cr to ^^Cr, the ex- 
citation energy of the first excited 2^ state decreases and 
i?4/2, the ratio of the excitation energy of the 4^ state to 
that of the 2^ state, increases. These data seem to indi- 
cate that a quantum phase transition from the spherical 
to deformed shapes takes place near TV = 40. The mi- 
croscopic origin of the enhanced quadrupole collectivity 
toward = 40 has been actively discussed from vari- 
ous theoretical approaches: the Hartree-Fock-Bogoliubov 
(HFB) mean-field calculations using the Skyrme force M 
or the Gogny force [1] , the spherical shell- model 0, flol |. 
and the projected deformed shell model [ill- These cal- 
culations have clarified the important role of the neutron 
(7g/2 and single-particle levels in the emergence of 
the quadrupole collectivity near A^ ^ 40. Although the 
spherical shell model calculations reproduce the experi- 
mental data rather well, the character of the quadrupole 
deformation, especially, the distinction between the equi- 
librium shape and shape fluctuations around it is not 
sufficiently clear. 

In this paper, we investigate the nature of the 
quadrupole collectivity in low-lying states of the neutron- 
rich Cr isotopes 58-64(^j, using an approach that treats 
the quadrupole deformations as dynamical variables. 
Thus, the distinction of the equilibrium shape and shape 
fluctuations is transparent. The deformation energy 
curve with respect to the axial quadrupole deformation 
was obtained in the Skyrme HFB mean-field calcula- 
tion 0, which shows that the quadrupole instability 
occurs around A^ = 38 — 42. However, the deformed 
minima are extremely shallow in these nuclei, suggest- 
ing a transitional character. In such transitional situa- 
tions, one naturally expects that large- amplitude shape 



fluctuations play an important role in determining the 
properties of low-lying excited states. Therefore, we take 
the five-dimcnsional (5D) quadrupole collective Hamilto- 
nian approach p^ . which is capable of describing the 
large-amplitude quadrupole shape fluctuations associ- 
ated with the quantum shape transition. It enables us 
to treat a variety of quadrupole deformation phenomena 
(vibrational, spherical-prolate transitional, rotational, 7- 
unstable, triaxial, oblate-prolate shape-coexistent situa- 
tions, etc.) on an equal footing. Dynamical variables of 
the 5D quadrupole collective Hamiltonian approach are 
the magnitude and triaxiality of quadrupole deformation 
(/3, 7) and the three Eulcr angles. To explicitly indicate 
the two vibrational and three rotational degrees of free- 
dom, we call it (2-|-3)D, hereafter. The 5D collective 
Hamiltonian is characterized by seven functions: the col- 
lective potential, three vibrational inertial functions (also 
called vibrational masses), and three rotational inertial 
functions. To evaluate the inertial functions, the Inglis- 
Belyaev (IB) cranking formula has been conventionally 
used. However, it is well known that the contribution 
of the time-odd components of the moving mean field is 
ignored in the IB cranking formula, which leads to the 
overestimation of excitation energies [13, il^ ■ 

The constrained Hartree-Fock-Bogoliubov plus 
local quasiparticle random-phase approximation 
(CHFB+LQRPA) method [H] is a method which 
can overcome the shortcoming of the IB cranking 
formula. It is derived on the basis of the adiabatic 
self-consistent collective coordinate (ASCC) method 
[isl [TtI [18| by assuming that there is a one-to-one 
mapping from a point on the collective submanifold 
embedded in the large-dimensional time-dependent 
HFB phase space to a point in the (/3,7) deformation 
space. In the CHFB+LQRPA method, the inertial 
functions are derived by transforming the local canonical 
coordinates determined by the LQRPA normal modes 
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to the (/3, 7) degrees of freedom. This method has been 
successfuUy applied to several phenomena: shape coex- 
istence/fluctuation in Se and Kr isotopes [ll . Il6l . 
development of triaxial deformation in ^^°Mo [20|, and 
shape fluctuations in neutron-rich Mg isotopes plj . 
Use of the Skyrme energy density functional in solving 
the CHFB+LQRPA equations has also been initiated 
for the axially symmetric quadrupole Hamiltonian 
[2^ . In this paper, we solve the LQRPA equations 
with use of the pairing-plus-quadrupole (P+Q) model 
including the quadrupole pairing interaction. For 
the collective Hamiltonian quantized according to the 
Pauli prescription, we solve the collective Schrodinger 
equation to obtain the excitation energies, vibrational 
wave functions, i?2-transition strengths and moments. 

This paper is organized as follows. We recapitulate the 
theoretical framework in Sec. [Ill In Sec. Illli we present 
results of calculation for 58-66(^j. discuss the nature 
of quadrupole collectivity in their low-lying states. To 
examine the role of the shape fluctuations toward the 
triaxial shape (7 degree of freedom), we also carry out 
calculations imposing the axial symmetry on the nuclear 
shape. We call it (1-|-2)D calculation indicating that the 
one vibrational degree of freedom in the /3 direction and 
the two Euler angles describing rotations perpendicular 
to the symmetry axis are taken into account. Comparing 
results of such a restricted calculation with those of the 
(2+3)D calculation, we shall see the importance of the 
triaxial shape fluctuations. We carry out a similar calcu- 
lation also for ^^Fe to compare with ^'^Cr. We then dis- 
cuss similarities and differences of the quadrupole shape 
transition near ^''Cr with iV = 40 and that near '^^Mg 
with N = 20. Conclusions are given in Sec. IIVI 



II. THEORETICAL FRAMEWORK 

In this section, we briefly summarize the framework of 
our collective Hamiltonian approach. See Ref. [l^ for 
details. 



A. 5D quadrupole collective Hamiltonian 

The 5D quadrupole collective Hamiltonian is given by 
Hcoii = T;ib + r,ot + l^(/3,7), (1) 

(2) 
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(3) 



k=l 



where Tvib and Trot represent the vibrational and ro- 
tational kinetic energies, while V the collective poten- 
tial energy. The moments of inertia are parametrized as 
7) = Ap^Dk{P, 7) sin^ 7fc with 7fc = 7 " 27rfc/3. 



The collective potential and inertial functions are de- 
termined with the CHFB+LQRPA method as explained 
in the next subsection. Once they are determined as 
functions of (/3,7), we quantize the collective Hamilto- 
nian according to the Pauli prescription. The collective 
Schrodinger equation for the quantized collective Hamil- 
tonian is given by 



{Tvib + Trot + V}-^^im{P, 7, ^) = E^i'^o.iMiP, 7, 



(4) 



where 
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Here, R[I3,^) and M^(/3,7) are defined as 

i?(/3,7) =A(/3,7)^2(/5,7)^3(/3,7), (7) 

(8) 

The collective wave function ^'q/m(/3, 7, r?) is specified 
by the total angular momentum /, its projection onto the 
z-axis of the laboratory frame M , and a distinguishing 
the states with the same / and Al . It can be written as 
a sum of products of the vibrational and rotational wave 
functions: 



<i>^iK{P,l){n\IMK), (9) 



where 



{n\iMK) 



21 + 1 



167:^(1 + Sko) 



(10) 



'^MK ^^'•^ Wigncr rotation matrix and K is the pro- 
jection of the angular momentum onto the z-axis in the 
body-fixed frame. The summation over K is taken from 
to / for even / and from 2 to / — 1 for odd /. 

The vibrational wave functions in the body-fixed 
frame, 7), are normalized as 



d/3d7|$„,(/3,7)nG(/3,7)|^ 



1, 



(11) 
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where 



which can be transformed to 



|ci>„K/3,7)P^ E \'^o.ik{P,i)\\ (12) 

/eleven 

and the volume element |G*(/3,7)| is given by 

|G(/3, 7)1 = 4/3«W^(/3, 7)i?(/3, 7) sin^ 37. (13) 

The symmetries and boundary conditions of the collec- 
tive Hamiltonian and wave functions are discussed in Ref. 

B. The CHFB+LQRPA method 

We determine the collective potential and inertial func- 
tions with the CHFB-hLQRPA method. We first solve 
the CHFB equation 

<5((/)(/3,7)I^chfb(/3,7)I0(/3,7)) =0, (14) 

i^CHFB = H-Y. ^^^^^'^^ - E ^''"'^2!' (15) 

T m 

with the constraints on the particle numbers and 
quadrupole deformation parameters: 

(0(/3,7)|iVM|'/'(A7)) =^^o"^ {r = n,p) (16) 

(m = 0,2). (17) 



(</>(/3,7)l^^™^l0(/3,7)) 



- n(+) 



Here, denotes Hermitian quadrupole operators, D20 
and {D22 +-D2-2)/2 for m = and 2, respectively. Then, 
we solve the LQRPA equations for vibration on top of the 
CHFB states obtained above, 

<5(0(/3,7)|[^chfb(/3,7),Q^(/3,7)] 

-iA(/3,7)l'/'(A7)) =0, (18) 

5(0(/3,7)|[ffcHFB(/3,7),-^'.(/3,7)] 
-C,(/3,7)Q'(/3,7)l'^(/3,7)) =0, (* = 1,2). (19) 

The infinitesimal generators, (5*(/3,7) and Pi{f3,^), are 
locally defined at every point of the (/3,7) deformation 
space. The quantity Ci{(3,j) is related to the eigen- 
frequency uJi{l3,j) of the local normal mode through 
^iiP^l) = Ci{(3,'-f). It is worth noting that these equa- 
tions are valid also for regions with negative curvature 
{Ci{f3,j) < 0) where uji{f3,^) takes an imaginary value. 

We assume that a 2D collective submanifold is embed- 
ded in a large-dimensional TDHFB configuration space 
and define local canonical coordinates ((71,(72) on it. We 
also assume that there is a one-to-one mapping between 
((71,(72) and (/3,7). By performing a similarity transfor- 
mation of ((71,(72), the vibrational kinetic energy can be 
written, without loss of generality, as 



I E 9^'' 

1=1,2 



Tvih — 2^^oo[i'20^]^ 



in terms of the time derivatives of ■ The vibrational 
masses with respect to are given by 



rS( + )j2 

(21) 



M„ 



'(A 7)= E 



The partial derivatives can be easily evaluated as 



(22) 



dD. 
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dD. 



(+) 



(23) 



9(7* dq'- 



= [DiV, jP^Wn)] , (24) 

without need of numerical derivatives. Through the def- 
inition of (/3,7) in terms of (£'20'' -^22^'')i ^^e vibrational 
inertial functions (D/?^, D^^, D^-,,) with respect to (/?,7) 
are easily obtained from (Mqo, M02, -^22) [3- To select 
two collective normal modes among the LQRPA modes 
obtained by solving Eqs. (flSl) and ([T9|) , we employ the 
minimal metric criterion [16|. 

The rotational moments of inertia are calculated by 
solving the LQRPA equations for rotation on each CHFB 
state: 

,5(0(/3,7)| [i?cHFB,^fc(/3,7)] - \{Jk)-'ik |(/>(/3,7)) =0, 

(25) 

(</)(/?, 7)1 [^fc(/3, 7), 4'] [(/-(Z?, 7)) = , 

(26) 

where Ik and 4'fe(/3,7) represent the angular momentum 
and the rotational angle operators with respect to the 
principal axes associated with the CHFB state |(^(/?, 7)). 
This is an extension of the Thouless-Valatin equation [2J| 
for the HFB equilibrium state to non-equilibrium CHFB 
states. We call J^k{P,j) determined by the above equa- 
tions 'LQRPA moments of inertia.' 

We solve the collective Schrodinger equation to ob- 
tain excitation energies and vibrational wave functions. 
Then, electric transition strengths and moments are read- 
ily calculated (see Ref. [l^ for details). 



C. Details of the numerical calculation 

(20) In this study, we adopt a version of the pairing-plus- 

quadrupole (P+Q) model [l2j including the quadrupole 
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pairing interaction as well as the monopole pairing in- 
teraction. We take two harmonic-oscillator shells with 
Nsh=3, 4 and Ngh^'^, 3 for neutrons and protons, re- 
spectively. The single-particle energies are determined 
from those obtained with the constrained Skyrmc-HFB 
calculations at the spherical shape using the HFBTHO 
code [IBl. In these Skyrme-HFB calculations, we employ 
the SkM* functional and the volume-type pairing with 
the pairing strength Vq = —180 McV fm""^. The pairing 
strength has been adjusted such that the calculated neu- 
tron pairing gaps at the HFB minima reproduce the ex- 
perimental gaps in ^^"^''^Cr. The single-particle energies 
are scaled with the effective mass of the SkM* functional 
m* /m = 0.79, because the P-l-Q model is designed to be 
used for single particles whose mass is the bare nucleon 
mass. 

The parameters of the P-f-Q model are determined 
as follows. For ^^Cr. the monopole pairing strengths 
Gq{t = n,p) and quadrupole interaction strength x ^-rc 
determined to approximately reproduce the HFB equi- 
librium deformation and the paring gaps at the spher- 
ical and HFB equilibrium shapes. For the other nu- 
clei 58,60,64,66Qj, g.^^^ 66pg^ .^g ^ssume the simple mass 
number dependence according to Baranger and Kumar 
[il: G5 - and x' = ^ A'^^^ {b denotes the 
oscillator- length parameter). We follow the Sakamoto- 
Kishimoto prescription to determine the strengths of the 
quadrupole pairing [26| . We omit the Fock term as in the 
conventional treatment of the P-l-Q model. 

The CHFB-I-LQRPA equations are solved at 60 x 60 
mesh points in the 7) plane defined by 

A = (i-0.5) X 0.01, (i = l,---60), (27) 
7, = (j-0.5)xl°, (.7 = 1, •••60). (28) 

For the calculation of the E2 transitions and moments, 
we use the effective charges (e|,g^ , e^^"* ) = (0.5,1.5). 



larger deformation in Cr, and the minimum becomes 
even deeper in ^''^Cr. In ^^Cr, the collective potential be- 
comes slightly softer than in ^^Cr. These potential energy 
surfaces indicate that a quantum shape transition from a 
spherical to a prolately deformed shape takes place along 
the isotopic chain toward = 40. In Fig. [U we plot the 
Nilsson diagrams of neutrons and protons as functions of 
(3 calculated for ^^Cr as in Ref. [l^- This is similar to 
Figs. 5(a) and 5(b) in Ref. Q. 

In Fig. [3l we plot the neutron and proton monopole 
pairing gaps Aq"^ (/3, 7) and Ag^^ (/3, 7), the vibrational in- 
ertial function £)^^(/3,7), and the rotational moment of 
inertia j7i(/3,7), calculated for ^^Cr. Figures[3] (c) clearly 
shows that the vibrational inertial function is well corre- 
lated with the magnitudes of the paring gaps: Dpp{j3,^) 

becomes small in the near spherical region where Aq"' 

and Ag^^ take large values. One might be concerned 
for complicated behaviors of D[)p{f3,^) in the strongly 
deformed region. However, they hardly affect low-lying 
states, because the collective potential energy is very high 
there and contributions from this region to the vibra- 
tional wave functions are negligibly small. Figures [3] (d) 
clearly indicates that the rotational moment of inertia 
also has a strong correlation with the pairing gaps. It 
takes the maximum value in the prolate region around 
P ~ 0.35. Both the neutron and proton pairing gaps 
become small there due to the deformed shell gaps for 
= 38 and Z = 24, see Fig. [2l In particular, the proton 
pairing gap almost vanishes. It results in the increase 
of the moment of inertia. As we shall see later, this en- 
hancement promotes the localization of the vibrational 
wave functions in the (/3, 7) plane for excited states with 
non-zero angular momenta. The rotational and vibra- 
tional inertial functions for the other isotopes are quali- 
tatively the same as those for ^^Cr. The enhancement of 
the moments of inertia mentioned above grows gradually 
with increasing neutron number up to A = 40. 



III. RESULTS AND DISCUSSION 

In this section, we present the numerical results for 
58-66Q]- discuss the nature of quadrupole collectiv- 
ity in their low-lying states. The results for ^^Cr are 
compared with those for the neighboring nucleus ^^Fe 
with A^ = 40. We furthermore discuss the similarities 
and differences with Mg isotopes around A^ = 20. 

A. Collective potentials and inertial functions 

We plot the collective potential V{j3, 7) calculated for 
58-66Qj. pjg [jj 'YliQ location of the absolute mini- 
mum is indicated by the (blue) circle. In ^^Cr, the ab- 
solute minimum is located at a nearly spherical shape. 
Although the minimum shifts to larger deformation in 
^"^Cr, the collective potential is extremely soft in the /3 
direction. A more pronounced local minimum appears at 



B. Yrast states in ss-es^j. 

We show in Fig. |3]the excitation energies of the 2^ and 
4^ states, their ratios i?4/2, the E2 transition strengths 
B{E2; 2^ 0^), and the spectroscopic quadrupole mo- 
ments of the 2^ states, together with the available ex- 
perimental data. The decrease in the excitation energies 
of the 2^ and 4+ states toward A' = 40 and the increase 
in their ratio from A^ = 36 to A^ = 40 are well described 
and indicate that the nature of the quadrupole collectiv- 
ity gradually changes from vibrational to rotational as 
the neutron number increases. However, the ratio R4/2 
at A^ = 40 is still 2.68, which is considerably smaller 
than the rigid-rotor value 3.33. The B{E2) values and 
spectroscopic quadrupole moments Q{2^) also suggest 
the onset of deformation: B{E2) increases and the mag- 
nitude of the spectroscopic quadrupole moments, which 
has a negative sign indicating a prolate shape, increase 
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with increasing neutron number and both of them reach 
a maximum at TV = 40. In Fig. HI we also plot the 
results of the (1+2)D calculations, in which only the axi- 
ally symmetric deformation is taken into account. Phys- 
ical meaning of the differences between the (2+3)D and 
(1+2)D calculations will be discussed in a subsequent 
subsection. 

In Table I, wc compare the results for ^"'Cr with those 
for ^^Fe. Experimental data indicate that the quadrupole 
collectivity is stronger in ^^Cr than in ^^Fe: the smaller 
E{2^) and the larger and B{E2) values for ^^-^Cr 

than those for ^^Fe. Our calculation reproduces these 
features quite well. 

We depict in Figs. [S] and [H] the squared vibrational 
wave functions multiplied by P'^ for the 0^,2^ and 
states in 58-66(^j. g^j^^^ those without the factor for the 
0^ and 2^ states in 6°Cr and ^''Cr, respectively. The 
factor comes from the volume element and carries its 
dominant (3 dependence (see Eqs. (fTT|) and (fT3|) ). The 
wave functions look quite different between the two cases. 
For instance, while the non-weighted 0^ wave function 
for ^'^Cr shown in Fig. [6] distributes around the spherical 
shape, the /3* factor changes it to the arcuate pattern seen 
in Fig. m In ^^Cr and ^OCr, the /J^-weighted 0;,^ wave 
functions exhibit arcuate distributions around (3 = 0.2 
covering the entire 7 region. Closely looking, one finds 
that, while the distribution for '^^Cr is almost uniform in 
the 7 direction, it is slightly leaning to the prolate side 
for 60Cr. 

With increasing neutron number, the 0^ wave function 
localizes more and more on the prolate side, reflecting the 
deepening of the prolate minima (see the collective poten- 
tial in Fig. [1]). In ^^Cr, the 0^ wave function still spreads 
over the entire 7 region, although it has a clear concen- 
tration on the prolate side. In ^^Cr, one can see a distinct 
peak around the prolate potential minimum, and the 0^ 
wave function is most localized at ^'^Cr. The vibrational 
wave functions clearly indicate the shape transition from 
spherical to prolate along the isotopic chain. 

For all these isotopes, one can see that the prolate peak 
grows with increasing angular momentum. This is due to 
the enhancement of the moments of inertia on the prolate 
side wc have already seen in Fig. [31 Even in ^^Cr whose 
ground state is rather spherical, the 2^ and 4^ states are 
weakly localized on the prolate side, which results in the 
finite spectroscopic quadrupole moment shown in Fig. |4l 



TABLE I: Excitation energies of the 2f state -E(2j^) in keV, 
the ratios _R4/2 of -5(4^) to -E(2j), and B{E2;2+ -> 0^) in 
Weisskopf units for ®'*Cr and ''^Fe. Experimental data are 
taken from Ref. [I,!!!]. 
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R4/2 


B{E2) 
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-R4/2 


B{E2) 


"''Ct 386 


2.68 


28.8 


420 


2.69 




®®Fe 685 


2.29 


15.5 


573 


2.47 


21.0 



In ^^Cr, the 2f and 4^ wave functions arc well localized 
on the prolate side, although the ground state wave func- 
tion still exhibit non-negligible shape fluctuation in the 7 
direction. Due to the growth of localization of the wave 
functions, higher angular momentum states acquire more 
rotor-like character than the ground state. This fact can 
be quantified by calculating the ratio 

i?6/4/2 ^ (EiGt) - E{2t))/{E{At) - E{2t)). (29) 

For instance, i?6/4/2 = 2.42 for ^"'Cr, which is fairly close 
to the rigid-rotor value 2.57, although the calculated Ri/2 
is 2.67 which is far from the rigid-rotor value 3.33. These 
results clearly indicate the importance of dynamical ef- 
fects of rotation on the nuclear shape. 

Lenzi et al. evaluated the intrinsic quadrupole 
moments Qint{I) for the yrast states of 62-66(^j- ^gjjjg 
the spectroscopic quadrupole moments Q{I) obtained in 
their shell-model calculation and the well-known relation 
between them for the axially symmetric deformation with 
K=Q. The resulting Q\nt{I) stay approximately constant 
along the yrast sequences in 62,64,66Qj.^ g^j^^j they inter- 
preted this as a fingerprint of a rigid rotor behavior. We 
have evaluated (5int(^) in the same way as Lenzi et al. 
but using our calculated Q{I). The resulting Qint{I) val- 
ues are similar to those of Lenzi et al. We feel, however, 
that this fact is insufficient to conclude that 62-66(^j, 
good rotors because the Q{I) are average values that are 
insensitive to the shape fluctuations. We need to exam- 
ine the properties of non-yrast states which are sensi- 
tive to shape fluctuation effects. We also note that Q{I) 
does not carry direct information about the ground state, 
and that, according to our calculation, the ground-state 
vibrational wave function is significantly different from 
those of the other yrast states with / ^ 0. 



C. Non-yrast states in ^»-^^>Cr 

To understand the nature of quadrupole collectivity, 
it is important to examine the properties of the excited 
bands including their interband transitions to the ground 
band, although they have not been observed experimen- 
tally yet. As typical examples of the calculated results, 
we display in Figs. [7]and[5]the excitation spectra and the 
B{E2) values of the low-lying states in ^°Cr and ^^Cr. 

Let us first discuss the ^°Cr case. We notice that 
the calculated excitation spectrum exhibits some features 
characteristic of the 5D harmonic oscillator (HO) limit; 
approximately equal level spacing in the ground band, 
approximate degeneracy of the 4^ and 2^ states, nearly 
equal values of B{E2;0^ 2f) and B{E2;4:t 2f), 
which arc about twice of B{E2;2f Of), etc. On the 
other hand, we also notice significant deviations from 
that limit. First, the Oj state is considerably lower than 
the 4^ and 2^ states. Second, the E2 transitions for- 
bidden in the HO limit arc sizable; e.g., those from the 
2^ state to the A'^ and 0^ states are fairly large. Third, 
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the B{E2) value from the 2^ state to the 2^ state is less 
than half of those from the 4^ and states. 

To examine the origin of these anharmonicities, let us 
look into the vibrational wave functions of the excited 
states. The /3^-weighted and non-weighted vibrational 
wave functions of the and 2^ states are displayed in 
Figs. [9] and [101 respectively. The wave function ex- 
hibits two components: one around the spherical shape 
and the other around (i = 0.35. Although it has a /3- 
vibrational feature, i.e., a node in the /3 direction, it 
also exhibits a considerable deviation from the 5D HO 
limit, in which the deformed component concentrating 
on the prolate side would spread uniformly over the 7 
direction. We can see a deviation from the 5D HO limit 
also in the 2^ state. The /?*-weighted 2^ wave function 
spreads from the prolate to the oblate sides. However, the 
non-weighted wave function reveals that it also has the 
/3- vibrational component. In fact, this state is a super- 
position of the large-amplitude 7-vibrational component 
spreading over the entire 7 region and the /3-vibrational 
component. In the 5D HO limit, the 2^ wave function 
has a node in the /3 direction, while the 2 J wave function 
has no node. The calculated 2^ wave functions indicate 
significant mixing of these components. Thus, the low- 
lying states in ^*^Cr may be characterized as a quadrupole 
vibrational spectrum with strong anharmonicity. 

Let us proceed to the ^^Cr case. We immediately no- 
tice some features different from ^°Cr. First, the ap- 
proximate degeneracy of the 4f and 23" states seen in 
^"Cr is completely lifted here. Second, the E2 transi- 
tions within the ground band are much stronger than 
those in ^°Cr. Third, two low-lying excited bands ap- 
pear: one consisting of the Q2 ,'^2 4^ states (ex- 
cited band I), and the other consisting of the 2^]", 3^, 4;^" 
states (excited band II, the 4j]' state not shown here is 
at 2.84 MeV). One might be tempted to interpret these 
excited bands in terms of the conventional concept of the 
/? and 7 bands built on a well-deformed prolate ground 
state, but, in fact, they are markedly different from them. 
First, there is a strong mixing of the /3- and 7-vibrational 
components, as seen from strong interband E2 transi- 
tions between the two excited bands. Second, the calcu- 
lated ratio of the excitation energies relative to i5(0j), 
{E{A+) - E{Q+))/{E{2+) - £;(0+)), is 2.51, which is far 
from the rigid-body value. Third, the iiT-mixing effects 
are strong, e.g., the A' = {K = 2) components of the 2^ 
(2;^) and 4^ (4;^) wave functions are at most 60%. To 
sum up, although the prolate deformation is apprecia- 
bly developed in the low-lying states of ^''Cr, the large- 
amplitude shape fluctuations play a dominant role and 
lead to the strong /3 — 7 coupling and significant inter- 
band E2 transitions. 

In Fig. [Til we plot the vibrational wave functions at 
7 = 0.5° and the probability density P(/3) of finding 
a shape with a specific value of /3 for the ground and 
excited 0+ states in ^°~^"'Cr. Note that the probabil- 
ity density vanishes at the spherical shape because of 
the factor in the volume element. It is seen that. 



while the ground-state wave function for ^"Cr distributes 
around the spherical shape, those for ^^Cr and ^''Cr ex- 
tend from the spherical to deformed regions with /3 ~ 0.4 
(see Fig. [TIJa)). Accordingly, the peak of the proba- 
bility distribution moves toward larger /3 in going from 
^°Cr to ^"^Cr (see Fig. fTTTb)). Concerning the excited 0+ 
states, their vibrational wave functions exhibit two peaks; 
a large peak at the spherical shape and a small peak at a 
prolate shape (see Fig. [TlTc)V In the probability distri- 
bution displayed in Fig. fTTT d). the spherical peaks moves 
to the /3 « 0.2 region and the peaks at /3 = 0.35 — 0.4 in 
turn become prominent. 

The above results indicate that large-amplitude shape 
fluctuations play an important role both in the ground 
and excited 0+ states. The growth of the shape fluc- 
tuations leads to an enhancement of the calculated EQ 
transition strengths p^{EO; 0^ — > O^,*') in going from ^^Cr 
to 62-66(^1.^ as displayed in Fig. [12l 

D. Comparison with the (1+2)D calculation 

To examine the role of dynamical shape fluctuations 
breaking the axial symmetry, we compare the above re- 
sult of the (2-|~3)D calculation with the (1-|-2)D calcula- 
tion where the axial symmetry is imposed. The proce- 
dure of the latter calculation is summarized in Appendix. 

The obtained excitation energies, E2 transitions and 
spectroscopic quadrupole moments are shown in Fig. |4l 
Obviously, the results of the (1+2)D calculation are more 
rotational than the (2+3)D results for all the isotopes; 
the i?4/2 ratios, the B{E2) values and the magnitude of 
Q{2'^) are larger while the excitation energies are smaller. 
The difference is most remarkable in the i?4/2 ratios. This 
implies that the neglect of the triaxial degree o f freedom 
leads to spectra much closer to the rotational limit. In 
fact, the (1-|-2)D calculation considerably overestimates 
the experimental i?4/2 values for ^''"^^Cr. The reason of 
this overestimation is easily understood in the following 
way. Let us first consider the /3-rigid limit where /3 is 
fixed at a certain non-zero value, (3 = /3o. In this limit, 
the (1-|-2)D calculation obviously gives the axial rigid- 
rotor value 3.33 for i?4/2. On the other hand, in the 
(2+3)D calculation, the shape fluctuation in the 7 degree 
of freedom is more or less always present. If the collective 
potential is flat in the 7 direction, as in the Wilets-Jean 
model 0, i?4/2 will be 2.50. In general, i?4/2 will take a 
value between 2.50 and 3.33 in the /3-rigid limit, depend- 
ing on the softness of the collective potential toward the 
7 direction (see Figs. 3 and 4 in Ref. [1^). In reality, 
as a quantum system, the nuclear shape always fluctu- 
ates in both the j3 and 7 directions. The experimental 
data that i?4/2 ~ 2.68 even in ^^Cr clearly indicate that 
such shape fluctuations cannot be ignored. In the ^°Cr 
case, the collective potential is quite soft in the /3 direc- 
tion and the (2-|-3)D calculation yields i?4/2 ~ 2.17. It is 
noteworthy that, even in such a situation, the 2^ and 4^ 
vibrational wave functions exhibit a localization about a 
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prolate shape (see Fig. [5]). This is due to the dynam- desired, 
ical effect minimizing the rotational energy for a given 
angular momentum. 
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E. Similarities and differences with the Mg 
isotopes around A*' = 20 

In Rcf. 01 , Lenzi et al. discussed similarities of the 
neutron-rich Cr isotopes near iV = 40 and the neutron- 
rich Mg isotopes around N — 20. Indeed, we have also 
found such similarities in our calculation. First of all, the 
growth of quadrupole collectivity in going from ^*'Cr to 
^''Cr is similar to that from ^^Mg to ^^Mg. In Fig. El 
while the ground state wave function in ^°Cr distributes 
around the spherical shape, they are considerably ex- 
tended to the prolately deformed region in ^*Cr. The 2f 
wave function has a peak on the prolate side in ^"^Cr and 
it shifts to larger [3 in ^^Cr. These features are similar 
to those seen in going from '^°Mg to '^^Mg in our calcula- 
tion [2l|. Concerning the excited Oj states in ^°Cr and 
^^Cr, as shown in Figs. fTIT c) and fTIT d). both vibrational 
wave functions exhibit a two-hump structure. Similar 
two-hump structures of the excited 0^ states have been 
obtained also in our calculation for "^^Mg and '^^Mg. 

On the other hand, we have also found significant dif- 
ferences between the ^''Cr region and the '^^Mg region. 
First of all, the K mixing is strong in the excited bands 
in the Cr isotopes, whereas it is weak in the Mg region. 
The shape fluctuations toward the 7 direction and the 
effect of the /3 — 7 coupling are larger in the Cr isotopes 
than in Mg. This can be clearly seen, for instance, in the 
2 J wave functions displayed in Figs. [9l andfTOl 



IV. CONCLUSIONS 

In this paper, we have investigated the nature of the 
quadrupole collectivity in the low-lying states of neutron- 
rich Cr isotopes ^^^^^Qr by solving the 5D collective 
Schrodinger equation. The vibrational and rotational in- 
ertial functions and the collective potential in the 5D 
quadrupole collective Hamiltonian are microscopically 
derived with use of the CHFB+LQRPA method. The 
calculated inertial functions include the contributions 
from the time-odd components of the moving mean field. 
The results of calculation are in good agreement with the 
available experimental data. The prolate deformation re- 
markably develops along the isotopic chain from A'^ = 36 
to 40. It is not appropriate, however, to characterize the 
low-lying state of Cr isotopes around ^"'Cr in terms of the 
prolate rigid-rotor model: the excitation spectra are still 
transitional and the large-amplitude shape fiuctuations 
dominate in their low-lying states. The calculated ex- 
cited bands exhibit strong couplings between the /? and 7 
vibrational degrees of freedom. For close examination of 
the nature of quadrupole collectivity in these nuclei, ex- 
perimental exploration of their excited bands is strongly 
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Appendix 

Here we summarize the procedure of the (l-f2)D calcu- 
lation to which the (2-|-3)D calculation reduces when the 
axial symmetry is imposed on the intrinsic shape. The 
dynamical variables in this reduced model are /? and two 
Euler angles {61,62), which describe the magnitude of 
the quadrupole deformation and the rotational motion 
perpendicular to the symmetry axis, respectively. The 
classical collective Hamiltonian of the (1-|-2)D model is 
given by 



ul + V{P), (30) 



k=1.2 



where uji are rotational frequencies related to the time 
derivative of the Euler angles. We quantize it accord- 
ing to the Pauli prescription, and then the collective 
Schrodinger equation reads 



vib 



k=1.2 



2^(/3) 



V ^o.im{P.6u62) = 

Ec.I^a.IM{P. 6^,62), (31) 



with 

^vib - 



92 



2Mp{P) 



1 dMp{(3) 1 dJiP) 



d_ 

w 

(32) 



The collective wave function can be written as a product 
of the vibrational and rotational wave functions: 

*„/M(/3,f^) = $a/(/3)(r!|/A./0). (33) 

The collective Schrodinger equation in the intrinsic frame 

is 

Tvib + ^-^^JfJ^ + *«/(/?) = E^i<i>^i{P), (34) 

and the ortho-normalization condition is given by 

d/33>:,(/3)3>„,7(/3)|G(/3)|i/2 =<5„„, (35) 
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with |G(/3)| = Mp{p)J\(i). 

For the (1+2)D calculation with the P+Q model pre- 
sented in subsection 3.3, we have used the same parame- 
ters as used in the (2-|-3)D calculation. In Ref. [22|, the 
(1+2)D collective Schrodinger equation ([34)1 was solved 
in the range —0.4 < /3 < 0.6 without respecting the 



boundary condition that the vibrational wave functions 
should satisfy at the spherical point (see Ref. [l^l)- In 
this paper, we have solved the collective Schrodinger 
equation in the range < /? < 0.6 respecting the bound- 
ary condition at /3 = 0. A detailed account on this point 
will be given in a future publication. 
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FIG. 1: (Color online) Collective potential energy surfaces 
V(/3,7) for ^^"^''Cr. The regions higher than 5 MeV (mea- 
sured from the HFB minima) are colored rosy-brown. 
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FIG. 2: Nilsson diagrams for neutrons (upper) and protons 
(lower) in ®^Cr as functions of /3, calculated as in Ref. 
The levels with the positive (negative) parity are plotted with 
solid (dotted) lines. 




FIG. 3: (Color online) Neutron and proton monopole pair- 
ing gaps, Aq"'(/3.7), Ag^'(/3.7), vibrational inertial function 
D,3^(/3, 7), and rotational moment of inertia J\{j3,-y), calcu- 
lated for ^^Cr. 
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FIG. 4: (Color online) (a) Excitation energies of the 2f states 
for 5*-^'^Cr. (b) Excitation energies of the states, (c) Ra- 
tios of E{4f) to -E(2+). (d) Reduced E2 transition probabil- 
ities B{E2;2j — >■ 0^) in Weisskopf units, (e) Spectroscopic 
quadrupole moments of the 2f states. Experimental data are 
taken from Refs. [2,[li,ll. 
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/3 fl ;3 fl 



FIG. 6: (Color online) Vibrational wave function squared, |^'a/x(/3, 7)1'^, for the O^*" and 2|j^ states in ^°Cr and ^*Ci. The 
contour lines are drawn at every twentieth part of the majcimum value. 
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FIG. 7: (Color online) Excitation energies and B{E2) values for ''"Cr in comparison with experimental data. Values on arrows 
indicate B(E2) in units of e'^fm*. Only B{E2) values larger than 1 Weisskopf unit are shown. Experimental data are taken 
from Ref. 0]. 
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FIG. 8: (Color online) Same as Fig. [Tlbut for ^''Cr. The experimental data are taken from Ref. 0|. 
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FIG. 9: (Color online) Squared vibrational wave functions multiplied by 13^, /^^ X^x \'^aiK{l3,'y)\'^ , for the and states in 
^°Cr and ^^Cr. 




FIG. 10: (Color online) Same as Fig. [6lbut for O2" and states. The contour lines are drawn at every twentieth part of the 
maximum value. 
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FIG. 11: (Color online) (a) Vibrational wave function 
squared |$q=i,i=o,a'=o(/3, 7 ~ 0.5°)^ of the ground states in 
60-64qj,^ (b) Probability densities integrated over 7, — 
J d7|$c=i,/=o,K=o(/?,7)l'|G(/3,7)r'''- (c) Same as (a) but for 
the 0^ states, (d) Same as (b) but for the states. 
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FIG. 12: (Color online) EO transition strengths p^{EO;0 
Ot) calculated for '^^-'"^Cr. 



